All stable state plots, with only the maximum tolerance strain

microxanox_version_here <- "0.9.1"
microxanox_version_path <- paste0("microxanox_v", microxanox_version_here)
if (packageVersion("microxanox") < package_version(microxanox_version_here)) {
  stop(paste0("microxanox version needs to be at least ", microxanox_version_here))
}
print(paste0("Version of microxanox: ", microxanox_version_here))

## [1] "Version of microxanox: 0.9.1"

wait_time <- "1e+06"
ss_9s <- readRDS(here::here("data",
                            microxanox_version_path, 
                            "temporal_method",
                            "event_definition_2",
                            paste0("ss_data_-1strains_waittime", wait_time, "_event_definition_2.RDS")))



max_CB_var_gmax <- max(ss_9s$CB_var_gmax_s)
max_SB_var_gmax <- max(ss_9s$SB_var_gmax_s)
max_PB_var_gmax <- max(ss_9s$PB_var_gmax_s)

ss_9s <- ss_9s %>%
  mutate(standvar_CB_gmax = CB_var_gmax_s / max_CB_var_gmax,
         standvar_SB_gmax = SB_var_gmax_s / max_SB_var_gmax,
         standvar_PB_gmax = PB_var_gmax_s / max_PB_var_gmax)

1 CB

this_dat <- ss_9s %>%
  filter(
    SB_var_gmax_s == 0,
    PB_var_gmax_s == 0
  )
this_dat <- this_dat[order(this_dat$CB_var_gmax_s), ]

for(i in 1:20)
  print(fig_state_vs_o2diff_sidebyside_alternative(this_dat[i,]))

#fig_state_vs_o2diff(only_SB_variable[13,])

2 SB

this_dat <- ss_9s %>%
  filter(
    CB_var_gmax_s == 0,
    PB_var_gmax_s == 0
  )
this_dat <- this_dat[order(this_dat$SB_var_gmax_s), ]

for(i in 1:20)
  print(fig_state_vs_o2diff_sidebyside_alternative(this_dat[i,]))

#fig_state_vs_o2diff(only_SB_variable[13,])

3 PB

this_dat <- ss_9s %>%
  filter(
    CB_var_gmax_s == 0,
    SB_var_gmax_s == 0
  )
this_dat <- this_dat[order(this_dat$PB_var_gmax_s), ]

for(i in 1:20)
  print(fig_state_vs_o2diff_sidebyside_alternative(this_dat[i,]))

#fig_state_vs_o2diff(only_SB_variable[13,])

4 CB-SB

this_dat <- ss_9s %>%
  filter(
    PB_var_gmax_s == 0,
    abs(standvar_CB_gmax - standvar_SB_gmax) < 0.00001
  )
this_dat <- this_dat[order(this_dat$CB_var_gmax_s), ]

for(i in 1:20)
  print(fig_state_vs_o2diff_sidebyside_alternative(this_dat[i,]))

#fig_state_vs_o2diff(only_SB_variable[13,])

5 CB-PB

this_dat <- ss_9s %>%
  filter(
    SB_var_gmax_s == 0,
    abs(standvar_CB_gmax - standvar_PB_gmax) < 0.00001
  )
this_dat <- this_dat[order(this_dat$CB_var_gmax_s), ]

for(i in 1:20)
  print(fig_state_vs_o2diff_sidebyside_alternative(this_dat[i,]))

#fig_state_vs_o2diff(only_SB_variable[13,])

6 SB-PB

this_dat <- ss_9s %>%
  filter(
    CB_var_gmax_s == 0,
    abs(standvar_SB_gmax - standvar_PB_gmax) < 0.00001
  )
this_dat <- this_dat[order(this_dat$SB_var_gmax_s), ]

for(i in 1:20)
  print(fig_state_vs_o2diff_sidebyside_alternative(this_dat[i,]))

#fig_state_vs_o2diff(only_SB_variable[13,])

7 CB-SB-PB

this_dat <- ss_9s %>%
  filter(
    abs(standvar_CB_gmax - standvar_SB_gmax) < 0.00001,
    abs(standvar_SB_gmax - standvar_PB_gmax) < 0.00001
  )
this_dat <- this_dat[order(this_dat$SB_var_gmax_s), ]

for(i in 1:20)
  print(fig_state_vs_o2diff_sidebyside_alternative(this_dat[i,]))

#fig_state_vs_o2diff(only_SB_variable[13,])